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, Abstract 

The correlator of vector and nonsinglet axial-vector currents in the external 
electromagnetic field is calculated within the instanton liquid model of QCD vac- 
uum. In general the correlator has two Lorentz structures: longitudinal wl and 
transversal wt with respect to axial-vector index. Within the instanton model the 
. saturation of the anomalous wl structure is demonstrated. It is known that in the 

chiral limit the transversal structure wt is free from perturbative corrections. In 
this limit within the instanton model we calculate the transversal invariant function 
wt at arbitrary momentum transfer q and show the absence of power corrections 
to this structure at large q 2 . Instead there arise the exponential corrections to wt 
at large q 2 reflecting nonlocal properties of QCD vacuum. The slope of wt at zero 
virtuality, the QCD vacuum magnetic susceptibility of the quark condensate and its 
momentum dependence are estimated. 
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1 Introduction 



Since discovery of anomalous properties PjJ |2] of the triangle diagram (Fig. 1) with 
^ . incoming two vector and one axial- vector currents [3] many new interesting results have 
been gained. Recently the interest to triangle diagram has been renewed due to the 
problem of accurate calculation of higher order hadronic contributions to muon anomalous 
magnetic moment via the light-by-light scattering process 1 . At low energies the dynamics 
of light-by-light scattering is nonperturbative, so one needs rather realistic QCD inspired 
model to find a solution with the lowest model sensitivity. 

The light-by-light scattering amplitude with one photon real and another photon has 
the momenta much smaller than the other two, can be analyzed using operator product 
expansion (OPE). In this special kinematics the amplitude is factorized into the ampli- 
tude depending on the largest photon momenta and the triangle amplitude involving the 
axial current A and two electromagnetic currents (one soft V and one virtual V). The 
corresponding triangle amplitude, which can be viewed as a mixing between the axial and 
vector currents in the external electromagnetic field, were considered recently in El • It 
can be expressed in terms of the two independent invariant functions, longitudinal wl and 

1 See, e.g., and references therein. 
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transversal wt with respect to axial current index. In perturbative theory for massless 
quarks (chiral limit) one has for space-like momenta q (q 2 > 0) 



w L (q 2 ) = 2w T (q 2 ) 
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(1) 




The appearance of the longitudinal structure is the consequence of the axial Adler-Bell- 
Jackiw anomaly PQ |2] . Because there are no perturbative (Fig. lb) [7j and nonperturba- 
tive (Fig. lc) [8 corrections to the axial anomaly, the invariant function wl remains intact 
when interaction with gluons is taken into account. Nonrenormalization of the longitudi- 
nal part follows from the 't Hooft consistency condition [S], i.e. the exact quark-hadron 
duality. In QCD this duality is realized as a correspondence between the infrared singu- 
larity of the quark triangle and the massless pion pole in terms of hadrons. It was shown 
in jHj (see also [5|) that in nonsinglet channel the transversal structure wt is also free 
from perturbative corrections. OPE analysis indicates that at large q the leading nonper- 
turbative power corrections to Wt can only appear starting with terms ~ 1/q 6 containing 
the matrix elements of the operators of dimension six ^U]. Thus, the transversal part 
of the triangle with a soft momentum in one of the vector currents has no perturbative 
corrections nevertheless it is modified nonperturbatively. 



Figure 1: Quark triangle diagram, (a); perturbative gluon, (6), and four-quark condensate, 
(c), corrections to it. 

In the present work we analyze in the framework of the instanton liquid model ^T] 
the nonperturbative properties of the triangle diagram in the kinematics specified above 
(see Section 2 for further details). The model is based on the representation of QCD 
vacuum as an ensemble of strong vacuum fluctuations of gluon field, instantons. They 
characterize nonlocal properties of QCD vacuum ^21EUE]- The interaction of light u, d 
quarks in the instanton vacuum can be described in terms of effective 't Hooft four-quark 
action with nonlocal kernel induced by quark zero modes in the instanton field (Section 
3). The gauged version of the model \W\ ITTj meets the symmetry properties with 
respect to external gauge fields (Section 4), and corresponding vertices satisfy the Ward- 
Takahashi identities. Below in Section 5 we demonstrate how the anomalous structure wl 
is saturated within the instanton liquid model. We also calculate the transversal invariant 
function Wt at arbitrary q and show that within the instanton model at large q 2 there 
are no power corrections to this structure. The nonperturbative corrections to wt at 
large q 2 have exponentially decreasing behavior related to the short distance properties of 
the instanton nonlocality in the QCD vacuum. We also estimate the slope of transversal 
invariant function at zero virtuality. 




(a) 



(b) 



(c) 
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When light quark current masses, rrif (/ = u,d) , are switched on, additional OPE 
structures appear, with the leading one being of dimension four ~ mfqa a pq. Its matrix 
element between vacuum and soft photon state is proportional to the quark condensate 
magnetic susceptibility introduced in Ref. [THj . Using the expansion of the triangle am- 
plitude in inverse powers of momentum transfer squared we will derive an expression for 
the magnetic susceptibility in the instanton model and find its momentum dependence 
(Section 6). 

2 The structure of VAV correlator 

We will employ a tensor decomposition of the VVA triangle graph amplitude suggested 
originally by Rosenberg [3] for the general kinematics of incoming momenta 

T^x(qi, q 2 ) = Ai$e P(a ,x + A 2 q p 2 e ppi/X + A 3 q^q P q^£ paflX (2) 
+ Mq 2 q P q 2 e p(7 pX + A 5 Qi q p q 2 £ P auX + A 6 q%q p q%e pauX , 

where q\ and q 2 are the vector field momenta with corresponding Lorentz indices \i and 
v. The coefficients Aj = Aj(qi,q 2 ),j = 1, ...6 are the Lorentz invariant amplitudes. The 
vector Ward identities provide a gauge invariant definition of the A\ and A 2 amplitudes 
in terms of finite amplitudes A^, k = 3, 6, 

Ai = (qiq 2 ) A 3 + q 2 A 4 , A 2 = ( qi q 2 ) A 6 + q\A h . (3) 

In the specific kinematics when one photon (q 2 = q) is virtual and another one (qi) 
represents the external electromagnetic field and can be regarded as a real photon with 
the vanishingly small momentum q± only two invariant functions survive in linear in small 
qi approximation [T^j. It is convenient to define longitudinal and transversal with respect 
to axial current index amplitudes [0] 

w L (q 2 ) = 4vr 2 2 4 (g 2 ) , w T (q 2 ) = 4tt 2 (l 4 (g 2 ) + A 6 (g 2 ) ) , (4) 

where tilted amplitudes are A(q 2 ) = A (gi = 0, q 2 = g) . In terms of w invariant functions 
the VAV amplitude becomes 

f puX (qi, q 2 ) = ^-j [wt (qhiEppvX ~ q 2 qiq 2 £ P ^x + q 2 qiq 2 £ P pau) - w L q\q{qle ppav \ . (5) 

Both structures are transversal with respect to vector current, q%T vX = 0. As for the axial 
current, the first structure is transversal with respect to g^ while the second is longitudinal 
and thus anomalous. 

The amplitude for the triangle diagrams (Fig.QJ can be written as a correlator of the 
axial current j x and two vector currents j u and j p 

Tpux = - J d 4 xdV^"^ (0| T{ Ju (x)j p (y)f x (0)}\0) , (6) 
where for light u and d quarks one has 

jp = qipQq, fx = qi\i5T 3 q , 
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the quark field q^ has color (i) and flavor (/) indices, the charge matrix is Q = ~ (| + r 3 ) 
and the tilted current is for the soft momentum photon vertex. 

In the local theory the one- loop result for the invariant functions Wt and wl is 2 

l-ioop _ 9 l-ioop _ 2N c f 1 daa(l - a) 

W L — Z W T — — — / — — 2 5 I ' 

3 Jo a(l - a)q 2 + mj 

where the factor N c /3 is due to color number and electric charge. The analytical result 
for the triangle diagram with finite quark masses has been obtained in [2H] by dispersion 
integral method. In the chiral limit, rrif = 0, one gets the result (0) (with additional 
factor N c /3). 

When nonperturbative contributions to the triangle amplitude (Fig. lc) are taken 
into account it was shown in PU| by using the OPE methods that at large Euclidean q 2 
the difference between the longitudinal and transversal parts, wlt = wl — 2wt, starts 
in the chiral limit from leading, ~ l/<? 6 , power behavior. The power terms are expected 
to contribute only into the transversal function Wt- Below we demonstrate that within 
the instanton liquid model in the chiral limit all allowed by OPE power corrections to wt 
cancel each other and only exponentially suppressed corrections remain. 



3 The instanton effective quark model 

To study nonperturbative effects in the triangle amplitude T^x at low and high momenta 
one can use the framework of the effective field model of QCD. In the low momenta 
domain the effect of the nonperturbative structure of QCD vacuum become dominant. 
Since invention of the QCD sum rule method based on the use of the standard OPE it 
is common to parameterize the nonperturbative properties of the QCD vacuum by using 
infinite towers of the vacuum expectation values of the quark-gluon operators. From 
this point of view the nonlocal properties of the QCD vacuum result from the partial 
resummation of the infinite series of power corrections, related to vacuum averages of 
quark-gluon operators with growing dimension, and may be conventionally described in 
terms of the nonlocal vacuum condensates ^21 Elj- This reconstruction leads effectively 
to nonlocal modifications of the propagators and effective vertices of the quark and gluon 
fields. 

The adequate model describing this general picture is the instanton liquid model of 
QCD vacuum describing nonperturbative nonlocal interactions in terms of the effective 
action [TT]. Spontaneous breaking the chiral symmetry and dynamical generation of a 
momentum-dependent quark mass are naturally explained within the instanton liquid 
model. The nonsinglet and singlet V and A current-current correlators, the vector Adler 
function have been calculated in fTf\ l%H 122] i n the framework of the effective chiral 
model with instanton-like nonlocal quark-quark interaction JH]. In the same model the 
pion transition form factor normalized by axial anomaly has been considered in j2H] for 
arbitrary photon virtualities. 

We start with the nonlocal chirally invariant action which describes the interaction of 

2 Here and below the small effects of isospin violation is neglected, considering m/ = m u — md- 
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soft quark fields [TH] 



S = j d 4 x q T {x) [i^Dp -m f ] q z (x) + X -G P j ' d*X j ' J] 
■ f(x n ) [Q(X - Xl ,X)T P Q(X, X + x 3 )Q(X - x 2 , X)T P Q(X, X + x 4 )] 



where _D M = — iV^ (x) — i-y^A^ (a?) and the matrix product (1 ® 1 + Z7sr a ® 275 r a ) pro- 
vides the spin-flavor structure of the interaction. In Eq. (jHj) g 7 = (w, d) denotes the flavor 
doublet field of dynamically generated quarks, Gp is the four-quark coupling constant, 
and r a are the Pauli isospin matrices. The separable nonlocal kernel of the interaction 
determined in terms of form factors f(x) is motivated by instanton model of QCD vacuum. 

In order to make the nonlocal action gauge-invariant with respect to external gauge 
fields V?(x) and A*(x), we define in (jHJ) the delocalized quark field, Q(x), by using the 
Schwinger gauge phase factor 

Q(x,y) = Pexp \i f dz„ [V«{z) + 75^ (z)] T a | qj(y), Q(x,y) = Q\z,y)<y», (9) 

where P is the operator of ordering along the integration path, with y denoting the 
position of the quark and x being an arbitrary reference point. The conserved vector and 
axial- vector currents have been derived earlier in ^SJ El 122] • 
The dressed quark propagator, S(p), is defined as 

S-^p) = ip-M(p 2 ), (10) 

with the momentum-dependent quark mass found as the solution of the gap equation 

M(p>) = m, + / ^J^ WTWW) - (U) 

The formal solution is expressed as ^2] 

M(p 2 ) =m f + (M q - m f )f 2 (p 2 ), (12) 

with constant M q = M(0) determined dynamically from Eq. (jllj) and the momentum 
dependent f(p) is the normalized four-dimensional Fourier transform of f(x) given in the 
coordinate representation. 

The nonlocal function f(p) describes the momentum distribution of quarks in the 
nonperturbative vacuum. Given nonlocality f(p) the light quark condensate in the chiral 
limit, M(p) = M q f 2 (p), is expressed as 

Its n-moment is proportional to the vacuum expectation value of the quark condensate 
with the covariant derivative squared D 2 to the nth power 

(0 \!D 2n a\ 0) = -N f M(p 2 ) 
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The nth moment of the quark condensate appears as a coefficient of Taylor expansion of 
the nonlocal quark condensate defined as [12] 



C{x) = (0 



q (0) Pexp 



A M (z) dzp 



q(x) 



0) (15) 



with gluon Schwinger phase factor inserted for gauge invariance and the integral is over 
the straight line path. Smoothness of C(x) near x 2 = leads to existence of the quark 
condensate moments in the l.h.s. of (JTljl for any n. In order to make the integral in the 
r.h.s. of ((H)) convergent the nonlocal function f(p) for large arguments must decrease 
faster than any inverse power of p 2 , e.g., like some exponential 3 

/ (p) ~ ex P ( — const -p a ) , a > as p 2 — ► oo. (16) 

Note, that the operators entering the matrix elements in (fT^j) and (|T3|) are constructed 
from the QCD quark and gluon fields. The r.h.s. of ()14j) is the value of the matrix elements 
of QCD defined operators calculated within the effective instanton model with dynamical 
quark fields. Within the instanton model in the zero mode approximation the function 
f(p) depends on the gauge. It is implied [121 Ej that the r.h.s. of (jTljl corresponds to 
calculations in the axial gauge for the quark effective field. It is selected among other 
gauges because in this gauge the covariant derivatives become ordinary ones: D — > d, and 
the exponential in ()15j) with straight line path is reduced to unit. In particular it means 
that one uses the quark zero modes in the instanton field given in the axial gauge when 
define the gauge dependent dynamical quark mass. The axial gauge at large momenta 
has exponentially decreasing behavior and all moments of the quark condensate exist. In 
principle, to calculate the gauge invariant matrix element corresponding to the of l.h.s. of 
(fT3)l it is possible to use the expression for the dynamical mass given in any gauge, but in 
that case the factor p 2n will be modified for more complicated weight function providing 
invariance of the answer 4 . 

Furthermore, the large distance asymptotics of the instanton solution is also modified 
by screening effects due to interaction of instanton field with surrounding physical vacuum 
[2H I2S] • To take into account these effects and make numerics simpler we shell use for 
the nonlocal function the Gaussian form 

f(p) = exp (-p 2 /A 2 ) , (17) 

where the parameter A characterizes the nonlocality size and it is proportional to the 
inverse average size of instanton in the QCD vacuum. 

The important property of the dynamical mass is that at low virtualities passing 
through quark its mass is close to constituent mass, while at large virtualities it goes to 
the current mass value. As we will see in Sect. 5 this property is crucial in obtaining 
the anomaly at large momentum transfer. The instanton liquid model can be viewed as 
an approximation of large- N c QCD where the only new interaction terms, retained after 
integration of the high frequency modes of the quark and gluon fields down to a nonlocality 

3 Very similar arguments lead the author of to the conclusion that finiteness of all transverse 
momenta moments of the quark distributions guarantees the exponential fall-of of the cross sections. 

4 If one would naively use the dynamical quark mass corresponding to popular singular gauge then one 
finds the problem with convergence of the integrals in (|14|l . because in this gauge there is only powerlike 
asymptotics of M (p) ~ p~ 6 at large p 2 . 
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scale A at which spontaneous chiral symmetry breaking occurs, are those which can be 
cast in the form of four-fermion operators (jSJ). The parameters of the model are then the 
nonlocality scale A and the four-fermion coupling constant Gp. 



4 Conserved vector and axial-vector currents 



The quark- ant iquark scattering matrix (Fig. 2) in pseudoscalar channel is found from the 
Bethe-Salpeter equation as 



T P (q 2 



Gi 



1 - GpJ PP (q 2 



with the polarization operator being 

d A k 



J PP (q 



f(k) f 2 (k + q)Tr [S(k) l5 S(k + q) l5 ] 



(19) 




Figure 2: Diagrammatic representation of the Bethe-Salpeter equation for the quark-quark 
scattering matrix, T, with nonlocal instanton kernel, I. 



The position of pion state is determined as the pole of the scattering matrix 



det(l -GpJ PP {q 2 )) I 



0. 



(20) 



J) T a 



The quark-pion vertex found from the residue of the scattering matrix is (k' = k + q) 

K(k,k') = 97Tqq z l5 f(k)f(k 

with the quark-pion coupling found from 

_ 2 dJa (q 2 



l-q 



dq 2 



(21) 



where m, is physical mass of the 7r-meson. The quark-pion coupling, g 2 q , and the pion 



decay constant, /„., are connected by the Goldberger-Treiman relation, g n = M q j f w , which 
is verified to be valid in the nonlocal model as requested by the chiral symmetry. 
The vector vertex following from the model (JSJ) is (Fig. Ek) 



r^k, k') = 1 , + {k + k'),M^\k ) k') 



(22) 



where M^(k, k') is the finite-difference derivative of the dynamical quark mass, q is the 
momentum corresponding to the current, and k (k') is the incoming (outgoing) momentum 
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Figure 3: Diagrammatic representation of the bare (a) and full (b) quark-current vertices. 
Diagram (c) shows separation of local (fat dot) and nonlocal parts of the full vertex. 

of the quark, k! = k + q. The finite-difference derivative of an arbitrary function F is 
defined as 

^M = m^>. (23, 

The full axial vertex corresponding to the conserved axial-vector current is obtained 
after resummation of quark-loop chain that results in appearance of term proportional to 
the pion propagator (THj (Fig. Eb) 

m f GpJp (q 2 ) 



r 5 Jk,k')= w + 2^f(k)f(k' 



Jap (0) - 



1-G P J PP (q 2 



(24) 



+ (k + k^^AP (0) 



(f(k')-f(k)) 



2 



where we have introduced the notations 



J (2tt) 

J AP (q 2 ) = 4:N c Nf J ^f M D { {J) ^ M{l + q)M{l) - (26) 
The axial- vector vertex has a pole at 

q 2 = -ml = m c (qq) / f 2 (27) 



where the Goldberger-Treiman relation and definition of the quark condensate have been 
used. The pole is related to the denominator 1 — GpJpp (q 2 ) in Eq. while q 2 in 

denominator is compensated by zero from square brackets in the limit q 2 — > 0. This 
compensation follows from expansion of J(q 2 ) functions near zero 

J PP {q 2 ) = G P l + m c (qq) M~ 2 (0) - q 2 g- 2 + O (g 4 ) , 
J AP (q 2 = 0) = M (0) , J P {q 2 = 0) = (qq) M" 1 (0) . 
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In the chiral limit mf — the second structure in square brackets in Eq. (J24j) disappears 
and the pole moves to zero. 

The parameters of the model are fixed in a way typical for effective low-energy quark 
models. One usually fits the pion decay constant, /„-, to its experimental value, which 
in the chiral limit reduces to 86 MeV [2E]. In the instanton model the constant, f^, is 
determined by 



f 

J 7T 



No 
4vr 2 



du u 



M 2 (u) - uM{u)M'{u) + u 2 M'(u) 
D 2 (u) 



{28) 



where here and below u = k 2 , primes mean derivatives with respect to u: M'(u) = 
dM(u)/du, etc., and 

D (k 2 ) = k 2 + M{kf. 
One gets the values of the model parameters [22J 

M q = 0.24 GeV, A P = 1.11 GeV, G P = 27 A GeV~ 2 . (29) 



5 VAV correlator within the instanton liquid model 

Our goal is to obtain the nondiagonal correlator of vector current and nonsinglet axial- 
vector current in the external electromagnetic field (VAV) by using the effective instanton- 
like model (JBJ). In this model the VAV correlator is defined by (Fig. 4a) 



d A k 
(M 3 



(30) 



■Tr [r„ (k + gi ,k)S(k + Ql ) r 5 x (k + qi ,k- q 2 ) S(k- q 2 ) T u (k, k - q 2 ) S (k)] , 

where the quark propagator, the vector and the axial- vector vertices are given by (JTUJ), (j22j) 
and (J2U), respectively. The structure of the vector vertices guarantees that the amplitude 
is transversal with respect to vector indices 



Tp,\{qi, q2)qi= T^xiqi, ^2)^2 = 



(31) 



and the Lorentz structure of the amplitude is given by ®. 

It is convenient to express Eq. (|30J) as a sum of the contribution where all vertices are 
local (Fig. 4b), and the rest contribution containing nonlocal parts of the vertices (Fig. 
3c). Further results in this section will concern the chiral limit. 

The contributions of diagram 4b to the invariant functions at space-like momentum 
transfer, q 2 = q 2 , are given by 



9q 2 



d A k 



7T 4 D\D- 



M+Hkq) 



% L) (q 2 ) = —At' if) ■ 



(32) 
(33) 



where the notations used here and below are 



(k+q) (k-q) 
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Figure 4: Diagrammatic representation of the triangle diagram in the instanton model 
with dressed quark lines and full quark-current vertices (a); and part of the diagram when 
all vertices are local one (b). 



D-t 



D(kl), 



M± = M(kl), f ± = f(kl). 



At large q 2 one has an expansion 



A[ L) (q 2 - oo) 



Ae 
67T 2 



Q 



i a (L) a {L) 



with coefficients given by 

J L ) - 



6 (4) 



du^M {2u + M 2 («)) 



D 2 (u) 



i {L) 

l (6) 



q 6 

2 r°° 

3 7o 



(34) 



du 



uM 2 (u) (u + 2M 2 (u)) 



(35) 

It is clear that the contribution ()32|) saturate the anomaly at large q 2 . The reason is that 
the leading asymptotics of (}32|) is given by the configuration where the large momentum is 
passing through all quark lines. Then the dynamical quark mass M(k) reduces to zero and 
the asymptotic limit of triangle diagram with dynamical quarks and local vertices coincides 
with the standard triangle amplitude with massless quarks and, thus, it is independent 
of the model. 

The contribution to the form factors when the nonlocal parts of the vector and axial- 
vector vertices are taken into account is given by 



A 



INL) 



d A k 1 



-M 2(1) (fc+,fc_) 2 



(kq) 
q 2 



M 4 

2 



M + - -M' + k 2 ± 



-(kq) 



One has for the leading terms of large q 2 asymptotics 



Af L) - oo) 



No 
Qtt 2 



(NL) 
l (4) 



+ 



(NL) 
l (6) 



+ O (q- S ) 



with coefficients given by 



a 



a 



(NL) 
(4) 

(NL) 
(6) 



du- 



uM (u) 
D 2 (u) 



(M (u) - uM' («)) 



(36) 



(37) 



(38) 



00 , u 3 M(u)M'(u) 
du- 



D 2 (u) 
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In sum of two contributions both power corrections with coefficients and are 
canceled. To prove cancellation for the a^) coefficient one needs to use integration by 
parts. 

Summing analytically the local (}3*2*j) and nonlocal ()36|) parts provides us with the result 
required by the axial anomaly 



w L (q 2 ) =Att 2 A 4 (q 2 ) 



2N C 1 
3 q 2 



(39) 



Fig. 5 illustrates how different contributions saturate the anomaly. Note, that at zero 
virtuality the saturation of anomaly follows from anomalous diagram of pion decay in 
two photons. This part is due to the triangle diagram involving nonlocal part of the 
axial vertex and local parts of the photon vertices. The result (|39|) is in agreement 
with the statement about absence of nonperturbative corrections to longitudinal invariant 
function following from the 't Hooft duality arguments. Earlier this consistency has also 
been demonstrated within the QCD sum rules (23 |2Hj and within dispersion method W. 
considering the lowest orders of expansion of the triangle diagram in condensates. 



For A^ L ^ (q 2 ) invariant function one gets 



.4 



(NL) 



wo = -a 



-~M' 
3 4 



N c r d A k 1 
6g2 J vr 4 D\D_ 

2k 2 , Mi - M_ 



(M+ + M_ 



M + - M (M+ - M. 



+ 



2Nc 
9q 2 



kl-kl 

d 4 k y'.W . M k\ 
vr 4 D 2 D_ k 2 -k 2 _ 



k 2 + 2 



(kq)' 



k 2 

6 (kq) K - 




(40) 



[M+ - M_ - 2M' + (kq)] . 



Then, let us consider the combination of invariant functions which show up nonper- 
turbative dynamics 



w L r(q 2 ) = w L (q 2 ) - 2w T (q 2 ) 



A 4 (q 2 ) + 2A 6 (q 2 ) 



(41) 



From (}3*3*j) we see that the contribution to wlt from the triangle diagram 4b with local 
vertices is absent. In sum of A^ (q 2 ) and Aq (q 2 ) a number of cancellations takes place 
and the final result is quite simple 



4N C f d 4 k y/Ml 
~3q~ 2 



n 2 D 2 D_ 




M, 



3 ± 



y]W+M (1) (A; + , k-) - 2 (kq) M+v / M (1 ' \k + , k 



(42) 



The behavior of wlt(q 2 ) is presented in Fig. 6. 

In the above expression the integrand is proportional to the product of nonlocal form 
factors / (k + ) f (k-) depending on quark momenta passing through different quark lines. 
Then, it becomes evident that the large q 2 asymptotics of the integral is governed by 
the asymptotics of the nonlocal form factor / (q) which is exponentially suppressed 
Thus, within the instanton model the distinction between longitudinal and transversal 
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Figure 5: Normalized wl invariant 
function constrained by ABJ anomaly 
from triangle diagram Fig. 4a (solid 
line) and different contributions to it: 
from local part, Fig. 4b, (dashed 
line) , and from the nonlocal part (dash- 
dotted line). 




Q(GeV) 



Figure 6: Normalized wlt invariant 
function versus Q predicted by the in- 
stanton model from the diagram Fig. 
4a. 



parts is exponentially suppressed at large q 2 and all allowed by OPE power corrections 
are canceled each other. Recently, it was proven that the relation 

w LT [m f = 0] = 0, (43) 

which holds at the one-loop level gets no perturbative corrections from gluon exchanges 
The instanton liquid model indicates that it may be possible that due to anomaly 
this relation is violated at large q 2 only exponentially. 

We also find numerical values of the slope of invariant function wlt at zero virtuality 

^f^l (/w) = -3.8 GeV" 2 , (44) 
d( l 2 q i=o 

and its width 

r f r \ _1 

,2 



X 2 LT = J uw LT (it) du ' (yf w LT (u) du^j = 0.54 GeV 2 . (45) 



6 Magnetic susceptibility of quark condensate 

In this section we consider the leading power corrections to wl and wt resulting from 
inclusion of current quark mass, m/, into consideration. An appearance of this kind of 
power corrections is already clear from perturbative expression (|7jl. In OPE the leading, 
by dimension, correction to the invariant functions wt,l{q 2 ) is 

A Wl = 2Au> t = -^^, (46) 

where Kf are the matrix element of dimension 3 operators 

Of = -iq f a^ 75 q f , (47) 
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between the soft photon and vacuum states. Proportionality to mj in (|46|) is in corre- 
spondence with chirality arguments. 

In perturbation theory the matrix element Kf of the chirality-flip operator Of is pro- 
portional to rrif. Nonperturbatively, however, due to spontaneous breaking of the chiral 
symmetry Kf does not vanish at rrif = 0. It is convenient to introduce the magnetic 
susceptibility Xm normalized by the quark condensate [IB] 



Kf = -4vr (qq) \m ■ 

This representation emphasizes that magnetic susceptibility for the nondiagonal vector- 
axial-vector vector correlator in the external electromagnetic field plays the similar role 
as the quark condensate for the diagonal correlators of vector and axial-vector currents. 

In the instanton model the VAV correlator is given by (}3T)j) with the quark propagator, 
the vector and the axial- vector vertices defined by (JTUJ), and (JUJ), with the quark 
mass, rrif, being included. Keeping in the calculation only linear in the current quark 
mass terms one finds at large q 2 the correction at twist 4 level (|47)|) for the contribution 
of diagram Fig. 4b 

m ,2 \ _ 12m f N c f, u 2 M (u) 



AA^ (q 2 - oo) = --^-^ 2 / du-^fl, (48) 
v ' g 4 3 ix l J D 6 [u) 

and from nonlocal part 
Alf L) {q 2 - oo) = [ du^pl [-u + 3M 2 («) - AuM («) M' {u)\ . (49) 

The leading asymptotics linear in current mass for the invariant function Aq is given by 
the relation 

Al 6 (q 2 -> oo) = -^Al 4 (q 2 oo) (50) 

which is in accordance with OPE (jlEj) . 

Then, summing up contributions (j4l?|l and ()49j) and comparing the result at large g 2 
with OPE one gets the magnetic susceptibility in the form 

1 N c f u{M(u)-uM'(u)) 
Xm M = J du ^ , (51) 

where the quark condensate is defined in (|T3J) . 

Alternatively, to get (|5T|) we may simply calculate the matrix element of O a ^ (|47j) 
between vacuum and one real photon state and use Eq. ()22j) for the quark-photon vertex. 
In this way it is easy to show that the result ()51|) stays unchanged when one includes the 
vector meson degrees of freedom. Indeed, in the extended model |2*2*1 13T] the vector vertex 
gets contribution from vector p and u mesons in the form 

(«, 

where T a is a flavor matrix, f v (p), Jy (q 2 ) , Gy are the nonlocal form factor, the polar- 
ization operator and four-quark coupling in the vector channel, correspondingly. Due to 
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conservation of the vector current one has By (q 2 = 0) = and thus there is no contribu- 
tion to the magnetic susceptibility. 

It is easy also to derive the momentum dependence of the magnetic susceptibility 



Xm (q) 



d A k 



(0\qq\0) J 4vr 4 



--k 2 ± M^(k + ,k. 



M+ - ^ (M, - M_ 



(1 + B V (q 2 ) flf_) - 



(53) 



presented in Fig. 7. At large q the integral in (}53j) is proportional to the quark condensate 
providing the asymptotic result 



Xm (q -> oo) 



(54) 




1 2 3 4 5 6 

Q(GeV) 



Figure 7: The momentum dependence of the magnetic susceptibility of quark condensate. 

Recently, Eq. (|51|) has been obtained by more complicated way in [32J. Also note 
that the instanton model does not support use of the pion dominance for estimate of the 
magnetic susceptibility as it was attempted in j^j. The reason is that the pion pole in the 
axial vertex (J23j) is accompanied by the exponentially suppressed residue Jp (q 2 ). Thus, 
it does not contribute to the twist 4 coefficient. 

Given model parameters one finds numerical values for the quark condensate and 
the magnetic susceptibility 

(0 \qq\ 0) (Minst) = - (214 GeV) 3 , X m (Minst) = 4.32 GeV" 2 , (55) 

where ni QSt is the normalization scale typical for instanton fluctuations. To leading- 
logarithmic accuracy scale dependence of these values is predicted by QCD as 

(0 \qq\ 0) (ji) = L-^' b (0 \qq\ 0) (/*) , Xm (fi) = L-^")l b Xm M , (56) 

where L = a s (/x) /a s (fio), b = (11N C — 2nj) /3, jg q = —3Cp is the anomalous dimension 
of the quark condensate, 70 = Cp is the anomalous dimension of the chiral-odd local 
operator of leading-twist, Cp = 4/3. We may fix the normalization scale of the model by 
comparing the value of the condensate with that found in QCD sum rule at some standard 
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normalization point: (0 \qq\ 0) (//q = 1 GeV) = — (240 GeV) 3 . Then one finds L = 2.17 
that corresponds to the normalization point /ii nst ~ 0.5 GeV, with the QCD constant for 

three flavors being Aq CD = 296 MeV. The rescaled magnetic susceptibility calculated 
in the model will be 

Xm (/imst = 1 GeV) = 2.73 GeV" 2 , (57) 

which is in rather good agreement with the latest numerical value of Xm obtained with the 
QCD sum rule fit [SSI: Xm (/^sr = 1 GeV) = (3.15 ±0.3) GeV -2 . A phenomenology of 
hard exclusive processes sensitive to the magnetic susceptibility Xm , see |3Sj, will possibly 
help to fix its value. 



7 Conclusions 

In the framework of the instanton liquid model we calculated for arbitrary momenta 
transfer the nondiagonal correlator of the vector and nonsinglet axial-vector currents in 
the background of a soft vector field. In this case we find that at large momenta the 
nonperturbative power corrections are absent in the chiral limit for the transversal part 
wt of triangle diagram. The transversal part is corrected only by exponentially small 
terms which reflects the nonlocal structure of QCD vacuum. Within the instanton model 
the saturation of the anomalous, longitudinal vol structure is demonstrated explicitly. 
Using the instanton liquid model we also derive an expression for the quark condensate 
magnetic susceptibility and its momentum dependence. 

The author is grateful to A. P. Bakulev, N. I. Kochelev, P. Kroll, S. V. Mikhailov, A. A. 
Pivovarov, O. V. Teryaev for helpful discussions on the subject of the present work. The 
author also thanks for partial support from the Russian Foundation for Basic Research 
projects nos. 03-02-17291, 04-02-16445 and the Heisenberg-Landau program. 
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